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The QED trace anomaly is calculated at one-loop level based on the loop regularization 
method which is realized in 4-dimensional spacetime and preserves gauge symmetry and 
. Poincare symmetry in spite ol the introduction of two mass scales, namely the ultraviolet 

^ ' (UV) cut-off Mc and infrared (IR) cut-off /ig. It is shown that the dilation Ward identity 

which relates the three-point diagrams with the vacuum polarization diagrams gets the 
standard form of trace anomaly through quantum corrections in taking the consistent limit 
Mc — oo and /is = which recovers the original integrals. This explicitly demonstrates that 
the loop regularization method is indeed a self-consistent regularization scheme which is 
applicable to the calculations not only for the chiral anomaly but also for the trace anomaly, 
at least at one-loop level. It is also seen that the consistency conditions which relates 
the tensor-type and scalar-type irreducible loop integrals (ILIs) are crucial for obtaining a 
consistent result. As a comparison, we also present the one-loop calculations by using the 
usual Pauli-Villars regularization and the dimensional regularization. 

PACS numbers: ll.lO.Gh, 11.30. Qc, 12.20.Ds. 



I. INTRODUCTION 



Symmetry and symmetry breaking play important roles in elementary particle physics. Accord- 
ing to the Noether's theorem, if a system is invariant under a continue global transformation, there 
is a conserved current corresponding to this transformation. Although the Noether's theorem is an 
exact one, it comes into exist at the classical level. Sometimes the consequences of the Noether's 
theorem are violated by quantum corrections and are related with anomaly. A well known ex- 
ample of anomaly is the chiral anomaly which is due to Adler, Bell and Jackiw In the 
evaluation of the anomaly, an unavoidable procedure is to regularize the loop integrals, therefore 
a gauge symmetry-preserving regularization is crucial for the calculation of quantum correction 
induced anomaly. Meanwhile, a calculation of quantum correction induced anomaly provides a 
useful laboratory for the check of the consistence of a new regularization prescription. 



Recently, a new regularization scheme named loop regularization was proposed 



It has 



been explicitly proved at one-loop level that the loop regularization can preserve the non-ableian 

n fl 

gauge symmetry [31 as well as supersymmetry l6[. Adopting such a symmetry-preserving loop 

II 

regularization, we have studied the chiral anomaly systematically [^. It was found that the two 
scales in the loop regularization, UV scale Mc and IR scale fis, play important roles in understanding 
the chiral anomaly. Under the condition Mc — )• oo which is applicable for the renormalizability of 
QED and a symmetrical treatment of all the three currents, the chiral anomaly in the massless 
QED can consistently be obtained in case of Hs = 0, while the massless QED could be free from 
anomaly if keeping /x^ ^ 0. For the massive QED, the chiral anomaly can be obtained with fig = 0, 
while if taking Hs ^ m, the chiral anomaly would also be absent. In Ref. [8(], based on the loop 
regularization method, the radiatively induced Lorentz and CPT violating Chern-Simons term 
in QED was calculated and a consistent result is obtained when simultaneously combiriing the 
evaluation for the chiral anomaly with ensuring the vector current conserved. In Refs. ^, [ij, the 
loop regularization method has been applied to study the gauge theory coupled to gravitation. 

In this paper, we will focus on another well-known anomaly, namely trace anomaly, which causes 
the violation of the dilation Ward identity due to the scaling invariance. The trace anomaly problem 
was firstly studied in Ref. 



in the framework of scalar field theory, and it was subsequently 



analyzed in various aspects 
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141. In this paper, we shall calculate the trace anomaly in the 
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framework of QED by using the symmetry-preserving loop regularization [3|, Since the standard 
form of the trace anomaly in QED is known, for example Ref. {3], our study can be regarded as 
a further independent check of the consistency of the loop regularization prescription. This comes 
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to our main purpose in the present paper. We shall show that with the symmetry-preserving loop 
regularization one can obtain the standard form of the trace anomaly when taking Mc — )• oo and 
/is = 0, which is required for recovering the original Feynman loop integrals and also applicable to 
QED as it is a renormalizable and IR divergence free theory. It is also seen that the consistency 
conditions relating the tensor type and the scalar type irreducible loop integrals, which have been 
verified in the loop regularization method, are crucial for the anomaly evaluation. From our explicit 
calculations, we arrive at the conclusion again that the loop regularization is indeed a consistent 
symmetry-preserving regularization. 

The paper is organized as follows: In Sec. m we discuss the dilation ward identity. Sec. IIIII 
presents our detailed calculations for the trace anomaly based on the loop regularization method. 
In Sec. IIVI we give, for a comparison, the calculations by using the Pauli-Villars and dimensional 
regulations. The last section is our conclusions. 

II. THE DILATION WARD IDENTITY 

The dilation transformation is a rescaling of the space-time coordinate 

^ x'^' = Xxi", (1) 

with A is a real constant. Under the dilation transformation, a field <I>(x) transforms as 

5d^{x) = i[B,^{x)] = {x-d + d)^{x), (2) 

where d is the scale dimension of the field ^> (x) and classically, it is the mass dimension of the field 
^(x), i.e., one for scalar field, three half for fermion field and others for composite operators. D is 
the generator of the dilation transformation. 

It was proved that [3] the current generated by the dilation transformation is related to 
the improved symmetric energy-momentum tensor 9^^ by 

D^{x) = x'^e^,, (3) 

then the trace of is the divergence of the dilation current 

9^;{x) = d^D'^ix). (4) 

By using the canonical commutative relations, one can show that the dilation charge D{xo) = 
J d^xD^{xi,XQ) can be taken as the generator of dilation transformation, that is, 

[D{xq),<^{x, xo)] = -i{x ■ d + d)^{x). (5) 



k + p 





FIG. 1: One loop diagrammatical representation of II^i^ and A^j^. 

With respect to this commutation relation, we can derive a simple dilation Ward identity (see 
Appendix R|) 



where 



d 

- P ■ -Q^)'nf^i'{p, -p) = A^^{p,-p), 



^^u{p,-p) = i J dW^-"(T*(J^(x)J,(0))), 
which can be diagrammatically illustrated in Fig. [T] 



(6) 

(7) 
(8) 



III. TRACE ANOMALY CALCULATION WITH LOOP REGULARIZATION 



The Ward identity ([6]) is a classical one and it is violated by quantum corrections [13j ]. In this 
section, we will investigate this violation using symmetry-preserving loop regularization [s, 4|. Our 
discussion is specified in the framework of QED with lagrangian 



From this lagrangian, the symmetric energy-momentum tensor can be easily got as 



(9) 



+^gf,uF^^Fxa-F^Fx,. 
Then, by using the equation of motion, we get 



(10) 



(11) 



For studying the quantum corrections to the dilation Ward identity ([6]) , one should consider the 
diagrams depicted in Fig. [TJ For the vacuum polarization diagram, after some algebra, one has 



111/ 



(-) 



-67, 



It — m {Iji + i/)) — m 



4e^ 



dx 



2l2^^{m) - h{m)g^^ + 2x(l - x){p^g^u - Pfj.Pu)Io{m) 



(12) 



where 



h^^y{'m) 
him) 
Io{m) 



kli ky 



(27r)4 (fe2 _M2)2^ 
d'^k 1 

(2^)4fc2_M2' 

d^A: 1 



(13) 



(27r)4 (fc2-M2)2' 

with M-^ = — x(l — x)p'^. The number subscript of / stands for the degree of divergence of 
the integral and specifically, "2" means the corresponding integrals are quadratically divergent, 
and "0" means the corresponding integrals are logarithmically divergent. After applying the loop 
regularizationy], we yield the consistency condition 



and the regularized vacuum polarization two-point function 



R 



.2/2 



Pt^Pf) / dxx{l - x)lQ{m). 
Jo 



The explicit forms of l2'{iTi) and lQ'{m) were given as [3] 



tR 
-'2 



tR 

-'o 



167r2 

i 



167r2 



ln^-7. + l + y2(^. 



(14) 

(15) 

(16) 
(17) 



with jj"^ = + JVP, -f^ =jE = 0.5772 ■ 



and 



yoix) = I da 







yi{x) 



a 
1 + x 



y2{x) = yo{x) - 2/1 (x), 



,0] 



(18) 



and, roughly speaking, they play 



where fig and Mc are introduced in the loop regularization 
the roles of UV and IR cutoff, respectively. Because the QED is a renormalizable and IR divergent 
free theory, we can safely set Mc to be infinity and /Ug to be zero, i.e., Mc — >■ oo and fig = 0. 

With the above analysis, after a simply algebra, the left hand side (l.h.s.) of ([6]) can be simply 
written as 



{2-p- ^)n^.(P, -P) = ^2ie\p^g^, - p^p,)p^ ^ dxx^l - xfl% 



{m). 



(19) 



To get the above relation, we have used the following expressions 

-M2 



tR 



167r2 2/x- 
y-2{x) = 1 -e"^, 



1 /- M 



and relation 



/^(m) = -4x(l - x^'^/^s 



(20) 



Now let us turn to the three-point function ([8]) which is illustrated in Fig. [TJ Including the 
crossing diagram, we get the expression 



^tiiy{P, -P) 



rTr 



(27r)4^^ V''^^'']^-m^"'{^^ + i>)-m""{^ + i))-m 

-64x^(1 - x)m^{p^py - g^uP^)I-2{m) 
-|-16a;(l — x)(l — 2x)m^5^jyp^/_2("i)|- 



-m- 



(21) 



After applying loop regularization with the consistency condition 



(22) 



we arrive at 



+16ie^m^c/^^/ j dxx{l - x){l - 2x)I%{m). 
Jo 



(23) 



From this expression, at the first sight, one may think that, for the triangle graph evaluated by 
loop regularization, the gauge invariance is not preserved due to the second term. In fact, this is 
not the case if one realizes the following general relation for the Feynman parameter integral 

'1 



dx{l-2x)f[x{l-x)] = 0, 







(24) 



which can be proved by changing the integration variable from a; to 1 — x, where (1 — 2x) = {l — x)—x 
is an odd function under such a change. We then arrive at a gauge invariant result 

-1 



^au{P:-p) = -2,2ie^m\p^pu - g^,uP^) I dxx{l - x)lf^2(.^)^ 



„2_2 



rR 







(25) 



which is similar to Eq. ()19p . 



Takeing the limit Mc — )• oo and using the exphcit form of I^2i have 
(2 - P • ^)n^^(p, -p) = -:^{PfMPu - g^lvP^) dxx^{l - x) 



From these two equations, we obtain their difference 

A,.(P, -p)-{2-p- ^)n,.(p, -P) = --,{p,Pu - 9,.P ) dxx{l - ^) ^2 + ^2_^(l_^)^2 - 



As QED is an IR divergence free theory, we safely take iXg = 0, which leads to the following simple 
relation 

d 

^^uiP^ -p) - (2 - P • ^)n^^(p, -p) = —^{Pf,Pu - g^^uP^) 



d 



-{2-p--^)Uj^,{p,-p,m = 0). (27) 



Thus the dilation Ward identity at one-loop order is given by 



-)n>,-,)-(2-p.- 



A^,(p, -p) = (2 - p • j-)n^AP, -p)-{2-p- j^)I^Up^ -P, m = 0) 



{2-p- j^)^UP^ -P) - ^(Pf^P- - 9,.uP')^ (28) 

which demonstrates that the dilation Ward identity ([6]) is violated by quantum corrections and the 
anomaly arises from the quantum effects. One can also see that our conclusion is exactly the same 
as that given in Ref. jisl - 

It is interesting to note that the anomalous term in Eq. (|28|) actually arises from the photon 
polarization H^u- This can be seen explicitly in the case of the massless QED. In the massless case, 
the classical ward identity becomes 

^i^Ap. -p) = {2-p- ^)^Ap^ -p. m = 0) = 0. (29) 
While one-loop contribution to the vacuum polarization gives 

(2-p— )n« (p,-p,m = 0) = 32ie2(pV-W)/_^ dxx" {I - xf , (30) 

which leads, in the limits Mc oo and fis = 0, to the following result 

i2-p--^)n^AP,-P,m = 0) = 32ie\p^g^,-pf,p,)p^ dxx\l - x)'| - 2[-p^x{l - x)] I 



e 



2 
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It is exactly the anomaly in the l.h.s. of (j28p . 

In an operator form, the anomalous Ward identity (|28|) can be written as 

= rn^^ + ^(^_F^^F,-^)^ (32) 

with as the fine-structure constant of QED. This one- loop level result agrees with the all order 
result given in Ref. [ij]. To all orders in perturbation theory, it has been proved that the trace 



anomaly can be expressed in the following form 



= (i + 5{aem))moi^i; + ^^^^N{-F^,Fn, (33) 

where N{^F^yF^^'^) is a subtracted form of the operator ^Zg^^F^j^F^^ which satisfies Eq. (2.3) in 
Ref. [14] after composed operator renormalization, and Z3 is the wave-function renormalization 
constant of photon. 5(aem) = — 1 is the QED anomalous mass dimension function, and 

/3(aem) = ''T^ ^dm QED j3 function in the on-shell scheme, with mo and m as the bare and 

physical masses of electron, respectively. Here we would like to mention that generally in the one 
coupling constant theory the (3 function starts to depend on the renormalization scheme from the 
three-loop level, and at present the QED /3 function in the on-shell scheme is obtained at four-loop 
order in Ref. [15] . Eq. (j33p means that the trace anomaly is proportional to /3-function to all 
orders. From the above one- loop order result of trace anomaly, we can obtain the well-known 



lowest order QED /3-function |l7t as /3(aem) = 



IV. TRACE ANOMALY WITH OTHER REGULARIZATION SCHEMES AT 

ONE-LOOP LEVEL 

In Pauli-Villars regularization, the regularized lagrangian of QED is 

n 

C'^ = Y1 C^{^^7^(i5/. - eAf,)xlji - rmiJiij}. (34) 

i=0 

It is enough to choose n = 2 to make the whole theory finite, and Cj and rrij can be specialized as 

Co = 1) Ci = I, C2 = —2, 

rriQ = m^, = + 2A^, 7712 = + A^, (35) 

with A as a constant representing the UV cutoff which should be set to infinity in the final result. 
Then the regularized A^^ is 

2 f' i i i 

A^^(p,-p) = -2 / 77r-uTr[e7Mir- — ^7^ /^ , 1 ^ v,/ , I ]• (36) 

~^ J (2^) It-mi {Jp+p) -rrn (J+i))-mi 



After Feynman parametrization, a direct one-loop calculation gives 

^Up^-p) = ^9^.uY,Cimj-^{p^.Pu-g^.up')Y.CiW?^ / drc-^^— (37) 

1=0 i=0 i \ JF 

By using the conditions (j35p . the first term vanishes. It should be mentioned that in Ref. [13] the 
mass terms is discarded directly without considering these condition because the coefficient Ci is 
omitted there. Then after the Feynman parameter integral, we obtain 

2 2 2 

^%{P,-P) = -^{P^Pu-9f.uP^)Y,a'^imjAi-l), (38) 

where Ai is given by 



A (n') = 2 - (/ - 4m,V)V^ 

'^■^ ' (p4 _ 4m2p2)l/2 p2 + (^4 _ 4^2^2)1/2 ' ^ ^ 

As required by the Pauli-Villar regularization, to get the meaningful result of the original theory we 
should take limit A — )• oo in the expression (j38p . So that in this limit case, all the heavy regulators 
give vanishing contributions, and the result can be written as 

A^^(p, -p) = -^{Pi^Pu - W^)^(m2Ao - 1). (40) 



For the one-loop vacuum polarization diagram H^y{p, —p), it has been calculated in Ref. 

The result can be read 

g2 , ^2 p2 

'^jiuiP, -P) = Y^iPfiPu - 5mi//)| In ^ - 6 dxx{l - x) ln[l - ^a;(l - x)]|. 



(41) 



from which we can easily obtain 

r 

—n) = -(n,,n,, — a,.,,n^) 

p2 6 



(2 - p • ^)n^,(p, -P) = -^{Pi^Pu - 9^.uP^) [^(m^^o " 1) " J ■ (42) 



2 



Combing ([38]) and (02]), we finally get 

A^.(p, -p) = (2 - p • ^)n^.(p, -P) - ^iPf^P'^ - P^9f..)- (43) 

This means that the classical Ward identity ([6]) is violated by an anomalous term. Compared 
with Eq. (|28p . we find that the anomalous term calculated with Pauli-Villars regularization is the 
same as that obtained with loop regularization. In the massless QED case, we can get the same 
conclusion as in loop regularization. 

For the calculation with dimensional regularization [l^ . one can draw the same conclusion. 
In the calculation with loop regularization, it has been seen that the two consistency conditions 
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(jl4p and (j22p are crucial to get the consistent conclusion. In dimensional regularization, these two 
conditions are satisfied, therefore the conclusions of the forms (jl5|) and (|25|) can directly be reached. 
The difference is that, the explicit forms of I^{m) and I^2i''^) dimensional regularization can 
be resulted from that in loop regularization by taking the limits Mc —> oo and Hs = 0. We give the 
one-loop results with dimensional regularization here 

n^i.(P> -P) = j^iPt^Pi^ - P 9fiiy)- - ^{PfiPu - P g,iu) J ax x{l - xjln — 

2e2 i-i ^2 

A/L(p, -p) = 2^P>^P'^ - P^dp^u) / dxx'^il- x)— (44) 

^ TT'^ Jo m'^ — x(l — xjp-^ 

From these equations, the calculation of anomaly is straightforward. We notice here that the 
anomalous term only depends on the finite parts of the above expressions, so that no matter what 
subtraction scheme be used we can always obtain the following result 

(2 - p • ^n%{p, -p) - -p) 



e f 
■^^iP/iPi^ -P^g^iu) dx 



m? — x(l — x)p^ 



e 



2 



= -^{PfiPu - P Qt^u)- (45) 

This subtraction scheme independent anomaly is the same as the one-loop results calculated with 
loop regularization and with Pauli-Villars regularization, and it demonstrates the self-consistency 
of applications of different methods. 



V. CONCLUSIONS 



In this paper, we have calculated the trace anomaly with loop regularization at one-loop level. 
The explicit demonstration has shown that the trace anomaly can be consistently obtained under 
the conditions ~^ and ^u^ = which are necessary for recovering original theory and are also 
applicable for our present calculations as QED is a renormalizable and IR divergence free theory. 
It has also been seen that the consistency conditions which are the direct deductions of non-abelian 
gauge symmetry are crucial for the anomaly evaluation. As the loop regularization preserves these 
conditions 0,01, it is then applicable for the studies on both symmetries and anomalies of gauge 
theories. 

We would like to mention that in low orders the anomalous dilation Ward identity (|28|) can be 
rewritten in the form of the classical one with the dimension d substituted by d' which is a sum of the 

n 

classical dimension and the anomalous dimension ^]. In the language of renormalization, this can 
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be understood as follows: As the vacuum polarization can be incorporated into the normalization 
constant of photon wave function Z3 , the renormalized photon wave function has the same form as 
the bare one except a scale dependent normalization constant Z3. Any regularization scheme will 
introduce at least one dimensional parameter into the bare theory, such as in loop regularization, 
A in Pauli-Villar regularization and in dimensional regularization. It is exactly this dimensional 
parameter which breaks the tree-level dilation Ward identity. On the other side, it is well-known 
that an exact scale invariant theory should be massless, so that the scale violation of a massless 
theory may be understood as the emergence of mass. In this sense, the radiative corrections can 
be considered as a possible way to generate mass [20|]. 
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In this appendix, we shall derive the classical Ward identity eq.®. Let us begin with eq.([8]) by 
considering the property of time order product, have 
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Appendix A: Derivation of Classical Ward Identity Eq. ([6]). 



^txv{P, -P) 



I 



d^xd'ye''P-y{T*{ei{x)J^{y)Um 




+ j d''xd^ye'P-ye{xQ)9{yQ-x^){J^{y)ei{x)MQ)) 

+ [ d^xdSe'p-yei-xo)eiyo){My)Mo)o^ix)) 



+ / d^xd''ye'P-ye{-xo)e{xo-yo){MO)di{x)J^{y)) 



+ / d^xd'ye'P-yei-xo)eiyo-xo)9i-yo){MO)My)0^xix)) 
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+ J d''xd'ye'P-ye{xo)e{-yo){e^^ix)MO)J^{y)). (Al) 
with respect to the relation eq.(j31) and the vanishing of the surface term, we then yield 

A^.(p,-p) = J d^xd^ye'P ye{xo)e{xo -yo)e{yo){dMx)J^{y)MO)) 
+ j d^xd^ye'P-ye{xo)e{yo-xo){J^{y)dMx)MO)) 
+ j d^xd''ye'P y9{-xo)9{yo){J^{y)MO)dMx)) 
+ I d^xd^ye'Pyei-xo)e{yo - xo)e{-yo){MO)J^{y)dMx)) 
+ J d^xdSe'P ye{-xo)e{xo-yo){MO)daDa{x)J^{y)) 

+ [ d^xd^ye'p-ye{xo)e{-yo){dMx)Mo)My)) 

d''ye^Py9iyo)myo),My)]MO)) 

dSe'pyeiyo){Mymo),,uo)]) 

d^ye^Pyei-yo){MO)[Diyo),My)]) 

d''ye'P ye{-yo)mO), J,(0)] J^(y)). (A2) 
which is simplified, by considering the commutation relation eq.([5]), to be 
A^.(p,-p) = i2d J dSe'P-ye{yo){T*{J^{y)J,m) 

+i j d^ye'J'-y9{yo){{yd)My)MO))+i J dV^-^^(-yo)( J.(0)(y • 9) J,,(y^A3) 
Noticing the following identities that 

i J d^ye'P-yeiyo){{yd)My)MO)) = -ip ■ ^ j dSe'^-ye{yo){J,{y)Jum 

-U I d''ye'P ye{yo){My)MO)), (A4) 

and 

i J dSe'P-ye{-yo){yd){MO)J,{y)) = -ip ■ ^ j dSe'^-ye{-yo){UQ)J,{y)) 

-4i j d'ye'Py9{-yo){MO)Jf.iy)), (A5) 
we then obtain the tree-level dilation Ward identity 

^,.{P,-P) = {2d-A-p-^}i I dSe^P-y{T*{J,{y)MO))) 
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= {2-p.^}U^-{p,-p), (A6) 
where Il^'^{p, —p) is defined in eq.d?]) and we have used the classical dimension of to be (i = 3. 



[1] J. S. Bell and R. Jackiw, Nuovo Cim. A 60 (1969) 47. 

[2] S. L. Adler, Phys. Rev. 177, 2426 (1969). 

[3] Y. L. Wu, Int. J. Mod. Phys. A 18, 5363 (2003) | arXiv:hep-th/0209021] . 

[4] Y. L. Wu, Mod. Phys. Lett. A 19, 2191 (2004) [arXiv:hep-th/0311082] . 

[5] J. W. Cui and Y. L. Wu, Int. J. Mod. Phys. A 23, 2861 (2008) l arXiv:0801.2T99l [hep-ph]]. 

[6] J. W. Cui, Y. Tang and Y. L. Wu, Phys. Rev. D 79, 125008 (2009) |arXiv:0812.0'892l [hep-ph]]. 

[7] Y. L. Ma and Y. L. Wu, Int. J. Mod. Phys. A 21, 6383 (2006) |arXiv:hep-ph/0509083] . 

[8] Y. L. Ma and Y. L. Wu, Phys. Lett. B 647, 427 (2007) |arXiv:hep-ph/0611199| . 

[9] Y. Tang and Y. L. Wu, Commun. Theor. Phys. 54, 1040 (2010) [arXiv:0807.033T ' [hep-ph]]. 

[10] Y. Tang and Y. L. Wu. larXiv:10T270626 [hep-ph]. 

[11] C. G. . Callan, Phys. Rev. D 2, 1541 (1970). 

[12] S. R. Coleman and R. Jackiw, Annals Phys. 67 (1971) 552. 

[13] M. S. Chanowitz and J. R. EUis, Phys. Rev. D 7, 2490 (1973). 

[14] S. L. Adler, J. C. Collins and A. Duncan, Phys. Rev. D 15, 1712 (1977). 

[15] D. J. Broadhurst, A. L. Kataev and O. V. Tarasov, Phys. Lett. B 298, 445 (1993) 
[ arXiv:hep-ph/9210255] . 

[16] G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 44, 189 (1972). G. 't Hooft, Nucl. Phys. B 61 (1973) 
455. 

[17] C. Itzykson and J. B. Zuber, "QUANTUM FIELD THEORY," New York, USA: McGraw-Hill (1980) 

[18] C. G. . Callan, S. R. Coleman and R. Jackiw, Annals Phys. 59 (1970) 42. 

[19] J. D. Bjorken and S. D. Dreh, "RELATIVISTIC QUANTUM MECHANIS," New York, USA: McGraw- 
Hill (1964) 

[20] S. R. Coleman and E. J. Weinberg, Phys. Rev. D 7, 1888 (1973). 



13 



